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THE PROMPT OPTICAL/NEAR-INFRARED FLARE OF GRB 050904: THE MOST LUMINOUS TRANSIENT 

EVER DETECTED 

D. A. Kann 1 , N. Masetti 2 , S. Klose 1 

The Astronomical Journal, in press 

ABSTRACT 

With a rcdshift of z = 6.295, GRB 050904 is the most distant gamma-ray burst ever discovered. It 
was an energetic event at all wavelengths and the afterglow was observed in detail in the near-infrared 
bands. We gathered all available optical and NIR afterglow photometry of this GRB to construct a 
composite NIR light curve spanning several decades in time and flux density. Transforming the NIR 
light curve into the optical, we find that the afterglow of GRB 050904 was more luminous at early 
times than any other GRB afterglow in the pre-Swift era, making it at these wavelengths the most 
luminous transient ever detected. Given the intrinsic properties of GRB 050904 and its afterglow, we 
discuss if this burst is markedly different from other GRBs at lower redshifts. 
Subject headings: Gamma rays: bursts — Gamma rays: individual: GRB 050904 



1. INTRODUCTION 

The detection of an extremely bright prompt optical 
flash accompanying GRB 990123 (Akerlof et al. 1999) 
highlighted the possibility of using GRB afterglows as 
backlighting sources to probe the high-redshift universe 
and the reinonization era (Lamb & Reichart 2000; Inoue 
et al. 2004; Totani et al. 2006). This promise was finally 
fulfilled with the discovery of GRB 050904 by the Swift 
satellite (Cusumano et al. 2006a, b), lying at a redshift 
of z = 6.295 (Kawai et al. 2006). Not only was a bright 
NIR afterglow discovered for this GRB (Haislip et al. 
2006; Tagliaferri et al. 2005), a prompt flash contempo- 
raneous to the 7-ray emission was seen by the robotic 
25cm TAROT telescope (Boer et al. 2006). 

Had the prompt optical flash of GRB 990123 been at 
a redshift of z = 1 and had it been unextinguished by 
any dust, it would have peaked at an apparent mag- 



nitude' 3 R = 7.6 
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(Kann et al. 2006, henceforth 



K06). This corresponds to an absolute i?-band magni- 
tude M R = -36.5 (2.3 • 10 16 L Qfl , assuming the V - R 
color for a G2V star from Ducati et al. 2001), making this 
the most luminous optical source ever detected at that 
time. In this paper, we compile all available data on the 
afterglow of GRB 050904 and construct a synthetic light 
curve for a hypothetical perfectly ionized universe show- 
ing no Lyman dropout up to the redshift of the GRB. We 
then derive light curve parameters and use the spectral 
energy distribution (SED) of the afterglow to extrapo- 
late the afterglow light curve into the R band assuming 
z = 1. We find that the prompt flash of GRB 050904 was 
the most luminous optical/NIR transient ever detected, 
even exceeding the _R-band luminosity of GRB 990123 at 
peak by more than one magnitude. 

2. THE COMPOSITE LIGHT CURVE OF THE GRB 
050904 AFTERGLOW 



Thiiringer Landessternwarte Tautenburg, Sternwarte 5, D— 
07778 Tautenburg, Germany 

2 INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica di 
Bologna, via Gobetti 101, 1-40129 Bologna, Italy 

3 In this paper, we use WMAP concordant cosmology with Hq = 
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2.1. Data mining and fitting methods 

The optical/NIR afterglow of GRB 050904 was ob- 
served by different telescopes over almost five decades in 
time and more than six decades in flux density. Due to 
its very high redshift, it was not detected at wavelengths 
shorter than the Ic band. We compiled all data from the 
following papers: Haislip et al. (2006), Tagliaferri et al. 

(2005) , Boer et al. (2006), Kawai et al. (2006), Gendre 
et al. (2006) and Berger et al. (2006). All data were cor- 
rected for Galactic extinction assuming Em-V) = 0.060 
(Schlegel et al. 1998), and then transformed into AB 
magnitudes (Oke & Gunn 1983). To transform the mag- 
nitudes into flux densities, we used the zero points given 
in Haislip et al. (2006) to transform their data and stan- 
dard zero points for all other data. Data in Kawai et al. 

(2006) , Berger et al. (2006) and z' data from Tagliaferri 
et al. (2005) are already AB magnitudes. 

The highest data density is available in the J band. 
Therefore, we started with the J band to construct a 
synthetic light curve by achromatically shifting the other 
bands to the J-band light curve. This method is justi- 
fied, as no systematic color evolution was found in the 
data set. For the normal forward shock evolution of the 
afterglow (once the cooling frequency v c has passed the 
optical/NIR bands) such achromacy is expected. While 
the early steep-to-shallow transition at t = 0.35 days 
may be chromatic, we only have J band data at this 
time anyway. In the following, we will always assume an 
achromatic evolution of the afterglow. 

We fitted the J-band data with a broken power-law (for 
more details on the fitting procedure, see Zeh et al. 2006, 
henceforth Z06). As Haislip et al. (2006) have noted, the 
early J-band data before 0.3 days are brighter than the 
extrapolation of later data to these times, indicating the 
steeper decline from what could be a reverse shock flash. 
Initially, we excluded these early data from our fit. We 
used the derived light curve parameters from the fit to 
the J-band light curve as a reference light curve and 
fitted the ii-band data with this reference curve. We 
thus found a color J — H. Assuming J — H = const., 
we shifted the ff-band data to the J zero point, merging 
the light curves. We then fitted the combined JH light 
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TABLE 1 

Results of the composite YJHK light curve fitting 



fit a x 2 d.o.f. m c h otQ a\ 02 tbl c *62 ™l d n 2 

0.35 ±0.14 -10 

2.45 ±0.18 ■•■ 2.63 ±0.37 •■• 1.82 ± 1.02 
-0.5 

2.4 ±0.4 ■•■ 2.6 ±1.0 



a The first fit uses data from 0.1 to 3 days after the burst, while the second fit uses only data after 0.3 days. H06 denotes the results from Haislip 
et al. (2006), and T05 those of Tagliaferri et al. (2005). 
^The apparent J-band magnitude at the respective break time is denoted as m c (cf. Z06, Eq. 1). 

c The break at which the forward shock afterglow begins to dominate over the prompt flash emission is £51, while t^2 denotes the jet break time. 
Break times arc given in units of days after the burst trigger. 
d ni and 712 arc the break smoothness parameters of the first and second breaks, respectively. 



1 

2 

H06 
T05 



59.2 
59.7 



17 
17 



18.80 ±0.47 
20.65 ±0.17 



1.39 ±0.13 



1.36 



+0.06 

0.07 



0.92 ± 0.05 
0.85 ± 0.08 

0.82+°;°! 

0.72+°^ 



curve again and shifted the iJ-band data points in small 
increments until \ 2 is minimized, thus deriving the true 
J — H color and the light curve parameters of the joint 
light curve which spans a longer time range (0.3 to 12 
days in the observer frame) and gives a better constraint 
on the post-jet break (after 2.6 days) decay slope a 2 (it 
is F u oc t~ a ). The X-band data (where we make no 
distinction between K, K s and K 1 ) and the F-band data 
were added in the same way. These additional shifts were 
not larger than 0.03 mag. 

We transformed the HST NICMOS F160W (AB) data 
point from Berger et al. (2006) to H (AB) by using the 
conv_AB value from the hyper Z package 4 (Bolzonella et 
al. 2000), m AB (H) = m AB (Fl60W) + 0.069. As the 
final HST data point is corrected for host contribution, 
we also corrected the other data points in the same way. 
From the synthetic host spectrum presented in Berger 
et al. (2006), we derived for the host galaxy a J-band 
flux density of 0.055 /iJy, which transforms into Jab = 
27.0. We assumed conservative errors of 0.3 mag. The 
effect of host subtraction is very small, at most 0.03 mag, 
since the last afterglow data point at 5 days is still 4 
magnitudes brighter than the host galaxy. 

2.2. Fitting the light curve 

The YJHK composite light curve consists of an early 
steep decay (Haislip et al. 2006), which goes over into 
a typical afterglow decay, followed by a further break, 
which is identified as a jet break (Tagliaferri et al. 2005). 
We fitted the two breaks separately, using the equation 
from Z06 (fixing the host flux to zero as the data have 
been corrected for host contribution). For the first fit, 
we excluded all data beyond 3 days from the fit. The 
break smoothness parameter n (cf. Z06) had to be fixed 
to n = — 10 (a negative value, as the break is from steep 
to shallow). The jet break was fitted with n left as a 
variable in the fit. For this fit, data earlier than 0.3 days 
were excluded. We thus derived three decay constants, 
which we label an, (the steep early decay), a\ and a 2 
(the typical pre- and post-jet break decay). We then 
transformed all data from AB back to Vega magnitudes 
using standard zero points. 

The z' data (B6er et al. 2006; Haislip et al. 2006; Tagli- 
aferri et al. 2005; Kawai et al. 2006) are already affected 
by Lya damping. While we used the composite reference 

4 http://webast.ast.obs-mip.fr/hyperz/ 




-0.5 c ' 1 

1E-3 0.01 0.1 1 10 100 

t (days after the burst in the observer frame) 



Fig. 1. — Composite z'YJHK light curve of the afterglow of 
GRB 050904, as it has been constructed with the procedure out- 
lined in §2.1 and 2.2. Inverted triangles show significant upper 
limits. The solid line represents the fit with the two broken power- 
laws, the dotted lines are the la confidence interval of the fit. The 
vertical dashed lines indicate the times of the two light curve breaks 
at 0.35 and 2.63 days, and the three decay slopes are labeled. The 
residuals Am represent observed minus fitted magnitudes. The 
residuals have been zoomed to show the scatter in the data. The 
afterglow is seen to peak at J Ri 10, almost three magnitudes above 
the extrapolation of the fit from 0.2 days, assuming a simple power- 
law decay. 



light curve to find the z' — J color, we did not implement 
the z' data in the derivation of the light curve parame- 
ters. We also did not implement the subluminous Z mea- 
surement from Haislip et al. (2006). Boer et al. (2006) 
give the TAROT measurements both as flux density at 
9500 A and as Ic magnitudes. We used the flux density 
and transformed it to z' magnitudes. Comparing these 
with the Ic magnitudes of Boer et al. (2006), we find 
Ic — z' = 2.49 mag. This is in very good agreement with 
the li — z' color of 2.48 mag derived from late afterglow 
data at 1.2 days, where I 2 is the VLT FORS2 /-band fil- 
ter (for details on the I\ and I2 filters, see Tagliaferri et 
al. 2005). The z' — J color derived at late times was then 
used to shift all z' data points to the J zero point (as 
always, we assume achromacy, i.e., a constant spectral 
slope). The result is a final composite light curve which 
includes all data from the early prompt emission to the 
late HST detection. In the following, when we speak of 
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Fig. 2. — The spectral energy distribution of the GRB 050904 
afterglow. The three dust model fits are shown. While MW (solid) 
and LMC (dashed) fits arc indistinguishable, the slight upturn of 
the SMC fit (dotted line) is seen. The Lyman absorption cutoff is 
not modeled, the gray fit curves arc extrapolations. Only upper 
limits (downward pointing triangles) have been been found in the 
Rc, r' and V bands. The flux density scale is arbitrarily chosen. 



the J-band light curve, we always refer to this composite 
light curve. 

2.3. Results of the light curve fitting 

The results of the light curve fits using the YJHK 
light curve are given in Table 1, and the complete com- 
posite light curve is shown in Fig. 1. As expected, a\ 
derived from the two different fits (cf. §2.2) agrees within 
errors. There is some scatter in the light curve data 
which leads to the high values of x 2 . This is either due 
to additional measurement uncertainties or small devia- 
tions from a power-law decay, as it has been found that 
the X-ray afterglow was extremely variable (Cusumano 
et al. 2006a, b; Watson et al. 2006). The values we de- 
rived are in full agreement with those of other authors 
(Table 1). These light curve parameters are typical for 
afterglows (see Z06 for a compilation of light curve pa- 
rameters of all pre-Swift afterglows). The light curve 
steepening Aa = a 2 — a.\ = 1.6 ± 0.2 is high but also 
not extraordinary (Z06, their Fig. 3). We note that this 
is one of the few light curves that allow the smoothness 
of the jet break to be fitted as a free parameter, and the 
result is in full agreement with the potential relationship 
between ct\ and n found by Z06 (their Fig. 8). 

For our fit, the extrapolation to very early times meets 
the earliest TAROT detection at 0.002 days, similar to 
what Boer et al. (2006) find. However, the following 
plateau phase around 0.001 days and the flare are much 
brighter than the extrapolated light curve (1 and almost 
3 magnitudes, respectively). This implies that the decay 
after the optical flare must be steep. Taking into ac- 
count the upper magnitude limit at 0.09 days (Gendre et 
al. 2006), we find a/; are > 3. Such a steep decay implies 
that the flare could be either due to reverse shock emis- 
sion or internal shocks from late central engine activity 
(Boer et al. 2006; Gendre et al. 2006; Wei et al. 2006; 
Zou et al. 2006). Since this is even steeper than cto, an 
additional break must exist, but no data are available 
during this time span. In particular, we find that the 
afterglow peaks 5 at J « 10, which is among the bright- 
est NIR afterglows ever detected and is a first hint at its 

5 The exposure time of the data point is f40 seconds (Boer et 



TABLE 2 

Results of the spectral energy distribution fitting. 



source a dust model /3 Ay c 



* 


none 


1.00 ±0.09 





* 


MW 


0.99 ±0.10 


0.02 ± 0.08 


* 


LMC 


0.92 ±0.35 


0.05 ± 0.20 


* 


SMC 


1.31 ± 1.20 


-0.10 ± 0.40 


H06 


none 


i 9C -+0.14 





T05 


none 


1.25 ±0.25 





T05 


none 


1.2 ±0.3 





P06 


none 


0.3 ±0.6 






a *: This work; H06: Haislip ot al. (2006); T05: Tagliaferri ct al. 
(2005) (two results from two different photo- 2 fitting procedures), P06: 
Price ct al. (2006). 



extreme luminosity. 

3. THE LUMINOSITY OF THE PROMPT FLASH 
3.1. Results of the spectral energy distribution fitting 

To construct the SED, we used the colors J — H, J — K 
and Y — J that were derived in §2.1 (Y, J, H and K data 
are not affected by Lyman a damping) . The SED is fitted 
as described in K06. First, we find an observed spectral 
slope 0q. Then the SED was fitted with an additional 
dust extinction curve, using dust models for the Milky 
Way, the Large and the Small Magellanic Cloud (Pei 
1992). This gives us the extinction in the V band in the 
host galaxy frame A^f st and an intrinsic spectral slope 
0, corrected for the host extinction. 

To extend the SED further into the optical bands (in 
the observer frame), we used reported Rc,r',V upper 
limits. By shifting the composite light curve to these 
limits, we produced additional Rc — J, r' — J and V — J 
colors which are lower limits on these colors. Instead 
of a flat or downward curved SED, we actually found a 
slight upward curvature, mainly due to the Y-band data 
(which have large errors) being slightly too bright (Fig. 
2, this can also be seen in Figs. 2 and 3 of Haislip et 
al. 2006). This results in a fit with SMC dust, which is 
typically favored for GRB host galaxies (K06, and ref- 
erences therein), finding a negative extinction, which is 
unphysical. The best result is actually obtained for Milky 
Way dust, even though this docs not imply that the host 
galaxy of GRB 050904 has dust resembling that of the 
Milky Way. While SMC dust is characterized by strong 
far-ultraviolet extinction, MW dust has a much flatter 
extinction curve at these wavelengths. Indeed, Maiolino 
et al. (2004) have found dust with flat FUV extinction in 
a quasar at almost the same redshift as GRB 050904, and 
Stratta et al. (2005) find that the optical and X-ray after- 
glow of GRB 020405 can be explained by dust with sim- 
ilar properties, although this GRB lies at a much closer 
redshift. The characteristic 2175 A bump of MW dust is 
not visible due to the very low extinction value we derive. 
We thus conclude that the extinction in the host galaxy 
is very low and the dust is dissimilar to SMC- type dust. 
In the following, we will use the values derived from the 
MW dust fit, deeming this as a conservative approach. 

The results of the SED fitting are given in Table 2. 

al. 2006). Possibly, the afterglow was even brighter on shorter 
timescales. 
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t (days after burst in the observer frame assuming z=1) 

Fig. 3. — Extrapolated light curve of the afterglow of GRB 050904 (thick line) assuming a fully ionized universe with no Lyman dropout, 
corrected for extinction by dust in our Galaxy and in the GRB host galaxy and shifted to z = 1, in comparison with the afterglows of the 
Golden Sample of K06. The light curve has been extrapolated into the Rq band as described in §3.4. The thin lines are the lc errors of 
the afterglow data, combined from measurement and extrapolation errors. The afterglow of GRB 050904 is more luminous than any other 
GRB afterglow at early times, making it the most luminous optical/NIR transient ever detected. At peak, it was about six magnitudes 
brighter than the early light curve of GRB 021004, which was the brightest afterglow at late times. GRB 021004 was also six magnitudes 
brighter than the extrapolation of the early afterglow of GRB 041006. At 0.05 to 0.8 days, the afterglow of GRB 050904 is among the 
most luminous afterglows, but later on, it is comparable to most other afterglows. The slanted dotted lines (with decay slope a = 1.2) are 
meant to guide the eye. 



Once again, our fits are in good agreement with those of 
other authors, except for Price et al. (2006), who ascribe 
the low value they derive to possible short term fluctua- 
tions. Haislip et al. (2006) and Tagliaferri et al. (2005) 
also note that while source frame dust extinction may be 
present, it must be low. The hydrogen column density 
along the line of sight in the host galaxy derived via op- 
tical spectroscopy is large (log N^f st = 21.3, Kawai et al. 
2006). Coupled with the very low extinction, this implies 
a dust-to-gas ratio that is even lower than for the SMC. 



The result is similar to GRB 990123, and in this case 
also, the prompt flash may have burned dust along the 
line of sight (cf. Fig. 6 in K06, and references therein). 

3.2. The circum-burst environment of GRB 050904 

Knowledge of the light curve parameters and the in- 
trinsic spectral slope allows us to derive conclusions on 
the nature of the circum-burst environment via closure 
relations (Price et al. 2002; Berger et al. 2002). We find 
that using u\ and (3, a model is strongly favored where 
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the cooling break frequency v c lies redward of the source 
frame ultraviolet after 0.4 days. This is substantiated by 
the spectral slope of the late X-ray afterglow, which is 
identical to that in the ultraviolet within errors. Wat- 
son ct al. (2006) find f3 x = 0.88 ± 0.04, Cusumano ct 
al. (2006a,b) find f3 x « 0.9 ± 0.1. Furthermore, Frail 
et al. (2006) used their radio detections of the afterglow 
of GRB 050904 to perform broadband modeling and ar- 
rived at the same conclusion. In this case, one cannot 
draw any conclusion on whether the GRB environment 
is of constant density (ISM) or is shaped by a progeni- 
tor wind. The observed very soft jet break (712, Table 1) 
may be indicative of a wind environment (Chevalier & 
Li 1999), whereas Frail et al. (2006) use an ISM model. 
Finally, Gendre et al. (2006) propose that the afterglow 
passes through a termination shock somewhere before 0.5 
days in the observer frame, switching from an early wind 
medium to a late ISM medium. Given the location of v c 
redward of the optical/NIR bands, an electron spectral 
index of p = 2/3 = 1.97 ± 0.29 is found, which is typical 
for GRB afterglows (e.g., Zhang & Meszaros 2004, and 
references therein) . The standard fireball model predicts 
p = a 2 (Zhang & Meszaros 2004), consistent with the a.2 
we find within 1.4er (Table 1). 

3.3. Derivation of the Rc-band light curve 

In the observer frame Rc band, the afterglow was un- 
observable due to Lyman damping. Typically, GRB af- 
terglows have the highest data density in the Rc band, 
and all afterglows that were examined in K06 were -Rc- 
band afterglows. Because of the synchrotron nature 
of the afterglow radiation (described by a power-law, 
F v oc with (3 being the spectral slope), and as- 

suming /3 =const., we can construct the Rc-h&nd light 
curve from the J-band light curve. The spectral slope 
we derived from the observer frame NIR data can be ex- 
trapolated into the observer frame optical, allowing us 
to find a Rc — J color without the damping influence of 
Lyman a. For this purpose, we used the observed spec- 
tral slope that has not yet been corrected for host galaxy 
extinction, (3q. Since the host galaxy extinction is very 
low, the spectral curvature is negligible and we can use a 
simple power-law for the extrapolation. This power-law 
has an error from the fitting process, which we added in 
quadrature to the measurement errors of the data points. 

3.4. The Rc-band luminosity of the prompt flash 

In order to place the afterglow of GRB 050904 in the 
context of the afterglow sample discussed in K06, we 
shifted the afterglow to a rcdshift z = 1. Using the 
derived intrinsic spectral slope (3, the source frame ex- 
tinction 4y s * (both from the MW dust fit) and the red- 
shift, we computed the magnitude shift denoted dRc 
in K06. This shift takes into account the extinction 
correction, the redshift difference and cosmological k- 
corrections (Sandage 1988) when moving the source from 
its original redshift to z = 1. The magnitude shift cre- 
ates a further error (from A/3 and A Ay) which we again 
added in quadrature to the errors of the data. We find 
dRc = — 5.05^°;Qg mag, which is higher than any value 
found in K06. Applying this correction to the i?c-band 
afterglow, we finally arrive at our main result: had GRB 
050904 lain at z = 1 and been unaffected by any line- 
of-sight extinction (including Lyman a damping), the 



optical flare detected by TAROT would have peaked at 
Rc = 6.48^0 28: almost visible to the naked eye. This 
is more than one magnitude brighter than the reverse 
shock peak of GRB 990123 (the flux density is higher 
by a factor of 2.8), making the afterglow of GRB 050904 
the most luminous optical transient ever detected. Boer 
et al. (2006), neglecting fc-correction, find a peak flux 
density of 1080 mJy (in the V band) for the afterglow 
of GRB 990123 and 1300 mJy (in the I c band) for the 
afterglow of GRB 050904 after scaling it to the redshift 
of GRB 990123, but do not remark upon this further ex- 
cept to call the two events comparable. For (3 = 0.99, 
the absolute -Rc-band magnitude is Mr = —37.6 ± 0.3 
(6.4 • 10 16 Lqr), making the afterglow in the Rc band 
alone between 60 and 600 times brighter than the intrin- 
sic bolometric luminosity of the gravitationally lensed 
hyperluminous BAL quasar APM 08279+5255, the most 
luminous persistent source known (Ibata et al. 1999). 

The complete composite Rc~ band light curve is shown 
in comparison with other afterglows in Fig. 3. It is im- 
mediately apparent that the optical afterglow of GRB 
050904 was much brighter than any other afterglow over 
almost two decades in time. At the peak of the prompt 
flash, it was six magnitudes brighter than GRB 021004, 
which is the brightest afterglow at late times, and an- 
other six magnitudes brighter than GRB 041006, one of 
the faintest GRBs in the sample. This span is much 
larger than at one day after the GRB, where the after- 
glow luminosities tend to cluster (K06; Liang & Zhang 
2006; Nardini et al. 2006). At later times, though, it 
resembles other afterglows. At one and four days after 
the burst (Fig. 3), we find R c = 18.13 and R c = 21.50, 
respectively, which transforms to absolute magnitudes 
M B = -24.51 ± 0.20 and M B = -21.14 ± 0.20, respec- 
tively (cf. K06, their Figs. 9 and 10). This ranks it 
among the brightest afterglows at one day, but it is not 
the brightest, as it has also been noted by Tagliaferri et 
al. (2005). At four days after the burst, it lies close to 
the center of the luminosity distribution. 

3.5. The spectral slope during the flare 

Since no multi-color data for the prompt emission 
phase are available, we shifted the early optical emis- 
sion achromatically, assuming no spectral changes be- 
tween this phase and the late afterglow phase (where no 
color evolution is evident). The assumption of achro- 
macy for the first 100 seconds after a GRB is surely an 
oversimplification. Strong spectral changes are expected 
within this time (Sari et al. 1998). However, does this 
imply that the true luminosity of the flare was actually 
lower than what we have derived? 

The spectral slope of the late optical/NIR afterglow, 
consistent with the cooling break lying at longer wave- 
lengths than the optical, is very steep, (3 « 1. During 
the prompt emission, the forward shock afterglow is ex- 
pected to be in the fast cooling phase (Wu et al. 2005, 
and references therein). If the medium surrounding the 
GRB progenitor is of constant density (ISM model), the 
optical band lies between the cooling frequency v c (at 
longer wavelengths) and the peak frequency v m , the ex- 
pected spectral slope is [3 = 0.5. In case even the cool- 
ing frequency has not yet passed the optical bands, the 
spectral slope in the optical actually rises, with (3 = —1/3 
(Zhang & Meszaros 2004, and references therein). A sim- 
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ilar situation would be seen if the early optical emission 
is dominated by the low energy tail of the prompt emis- 
sion, a situation which was seen for GRB 041219A and 
GRB 050820A (Vestrand et al. 2005, 2006). This has 
also been suggested for GRB 050904 (Wei et al. 2006), 
although the optical flux in the flare is much brighter 
than the extrapolation of the X-ray spectrum for GRB 
050904 (Boer et al. 2006), in contrast to the cases of GRB 
041219A and GRB 050820A (Vestrand et al. 2005, 2006). 
We thus conclude that the spectrum in the optical bands 
during the prompt emission was not steeper than at late 
times, and possibly shallower, i.e. (3 < 1. 

On the one hand, this implies that the Rc magnitude 
of the afterglow at z — 6.295 would be brighter than what 
we derived with the spectral slope of the late afterglow. 
On the other hand, when we apply the cosmological shift 
to z = 1, the dRc shift would become smaller. We find 
dRc = —4.26 mag for j3 — 0.5 and dRc = —3.10 mag for 
[3 = —1/3 (also correcting for the small host extinction), 
resulting in Rc = 6.93 and Rc — 7.52, respectively. 
Thus, for j3 = —1/3, the extrapolated i?c-band magni- 
tude at z = 1 of the flare would be comparable to that 
of GRB 990123. However, we should also note that the 
same argument may apply to the flare of GRB 990123, 
so that the flare of GRB 050904 truly is brighter with a 
high probability. 

4. IS GRB 050904 DIFFERENT FROM OTHER LONG 
GRBS? 

In the literature on GRB 050904, consensus is not 
reached on whether GRB 050904 is, apart from its red- 
shift, a special event. Let us review some aspects of both 
the prompt emission and the afterglow: 

In many ways, GRB 050904 does not differ from the 
typical mean values found for GRBs. In the source frame, 
T 90 = 31 seconds, with 7-ray emission detected up to 69 
seconds after the trigger, including the flare seen promi- 
nently in the X-rays and the optical (Boer et al. 2006; 
Cusumano et al. 2006a, b). The intrinsic luminosities of 
the X-ray (Watson et al. 2006), the late optical (as shown 
in this work) and the radio afterglow (Frail et al. 2006) 
are typical, as are the light curve decay parameters and 
the jet break time in the rest frame (Z06, K06). 

But some aspects of GRB 050904 do stand out. The 
peak energy E p lies beyond the Swift BAT energy range. 
From joint fitting of Swift BAT and Suzaku WAM data, 
Sugita et al. (2006) report a preliminary E p = 366I/143 

keV, translating into E p = 26701^43 keV in the rest 
frame. This peak energy is the highest found in the com- 
pilation of 56 bursts (Amati 2006), only GRB 050717 
(observer frame E p — 2101 tg^ 4 keV, time-integrated 
spectrum, Krimm et al. 2006) very probably has a 
higher peak energy (the rcdshift is not known but as- 
sumed to be high). The isotropic energy is not ex- 
actly known, Cusumano et al. (2006a, b) derive Ei so = 
6.6 x 10 53 - 3.2 x 10 54 erg. The upper end of the 
range is higher than for any other burst, including GRB 
990123 (Amati 2006), and even the lower value is higher 
than for most GRBs with well-determined £i S o. Us- 
ing a jet break time of U = 2.6 ± 1.0 days, Tagliaferri 
et al. (2005) derive log £ 7 (crg)= 51.6 - 52.1. This is 
the second highest beaming-corrected energy ever found 
(Friedman & Bloom 2005), after GRB 050820A (log 



£; 7 (erg)= 51.7-52.2, Cenko et al. 2006).The X-ray after- 
glow is the most variable ever seen by Swift (Cusumano 
et al. 2006a,b; Watson et al. 2006), implying that, in 
the source frame, the central engine was active for sev- 
eral hours after the burst (Zou et al. 2006). We have 
shown in this paper that the early optical/NIR emission 
is brighter than for any other GRB detected so far. Frail 
et al. (2006) derive a high circum-burst density of 680 
cm~ 3 from broadband afterglow modeling, noting that a 
density typical for GRB environments (10 cm -3 ) is ruled 
out by the absence of a bright early radio afterglow. 

There is no aspect of GRB 050904 that exceeds the 
typical GRB parameters so much that one could claim 
that the progenitor of GRB 050904 was markedly differ- 
ent from those of all other long GRBs, e.g. a Population 
III star. However, it is clear that GRB 050904 is one of 
the most extreme GRBs ever observed, comparable only 
to GRB 990123 (cf. Boer et al. 2006), and exceeding even 
this seminal event in some aspects. 

5. CONCLUSIONS AND OUTLOOK 

We have shown that the early flare of GRB 050904 
detected by the TAROT telescope is the most luminous 
optical/NIR transient ever detected, exceeding even the 
peak magnitude of the prompt optical flash of GRB 
990123. While GRB 050904 is not markedly different 
from other GRBs in any single aspect, it is among the 
most extreme GRBs ever detected. 

So far, only a single GRB at the end of the reionization 
era has been detected. It is thus not possible to decide 
whether GRB 050904 is typical for GRBs at a very high 
rcdshift, or if it is simple observational bias that the first 
GRB to be found at z > 6 is one of the most extreme 
events ever detected. More typical events would not have 
sufficient follow-up to determine the necessary parame- 
ters, especially the redshift. The Swift satellite, with 
its high sensitivity in an energy band that is typical for 
highly rcdshifted hard bursts and its possibility of long- 
duration image triggers, is a very powerful instrument for 
detecting these very distant events. Probably, more very 
high redshift GRBs will be detected during the lifetime 
of Swift, possibly even breaking the redshift records for 
quasars (Fan et al. 2003) and galaxies (lye et al. 2006). 
Indeed, a GRB at a higher redshift than GRB 050904 
may have already been found, GRB 060116 (Grazian et 
al. 2006) 6 . 

While the detection of a GRB by Swift is the first and 
most important step, most parameters characterizing the 
event can only be obtained by ground-based telescopes. 
The determination of the true burst energetics depends 
on the measurement of the redshift and the jet break 
time, and GRBs at z > 6 are bright only in the NIR. 
While optical spectroscopy is still feasible to z ~ 7, any 
higher rcdshifts will move the metal lines that are used 
for an exact determination of z beyond the typical limit 
of optical spectroscopy (and the region beyond 0.7 /im is 
already strongly affected by sky lines). NIR spectroscopy 
even at large facilities needs bright sources, and thus far, 
NIR spectroscopy of a GRB afterglow has never been 
successful in detecting any lines. The dilemma is that 
observing time at large facilities is precious and observa- 

6 However, see also Piranomonte et al. (2006) and Tanvir et al. 
(2006). 
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tions would only be triggered in case there is already a 
strong suspicion that the event lies at a very high red- 
shift. This implies the need for rapid photometric red- 
shift determinations that use the Lyman a trough, as 
it has been successfully applied to GRB 050904 (Hais- 
lip et al. 2006; Tagliaferri et al. 2005; Price et al. 2006). 
One telescope that is already in use and has the capa- 
bility for rapid multi-color photometry is the MAGNUM 
telescope, which also observed GRB 050904 (Price et al. 
2006). A key element would be the deployment of more 
moderately large robotic NIR telescopes such as PAIRI- 
TEL (Bloom et al. 2005) and REM (Zerbi et al. 2001), 
which also has an optical camera. 
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More rapid follow-up of high rcdshift bursts will finally 
allow GRBs to fulfill their promise as the ultimate probes 
of the very early universe. 
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